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ДИНАМИЧЕСКИЕ МОДЕЛИ УПРАВЛЯЮЩИХ НАГРУЗОК 
РЕОЛОГИЧЕСКИХ СИСТЕМ 


Предложена модель поведения реологического тела (вязкоупругого материала, конструкции или систе- 
мы) при управляющем воздействии нагрузки, действующей по заданному закону в течение некоторого 
времени. 

Ключевые слова: модель реологической среды, наследственные уравнения, полимерные композитные 
материалы. 


Введение. Для моделирования поведения сложных динамических систем необходимо вводить 
соответствующие сложные функции безразмерных параметров, отражающих весь комплекс 
действующих на объект нагрузок. Чтобы получить такие зависимости, целесообразно исполь- 
зовать методы масштабного моделирования, подробно описанные в литературе [1]. 
Постановка задачи. В качестве управляющей нагрузки удобно рассмотреть функцию гармо- 
нического изменения действующих напряжений о(1), в качестве управляемой функции отклика 
системы -— соответствующую ей функцию деформаций, также изменяющуюся по гармониче- 
скому закону с запаздыванием по фазе, равным углу механических потерь Ло. 
Моделирование поведения реологической среды. Универсальной моделью поведения 
реологической среды являются наследственные уравнения типа Больцмана и Вольтерра для 
полимерных композитных материалов и подобных им реологических сред. Однако наиболее 
наглядно поведение систем при данных условиях можно описать с помощью динамических 
моделей вязкоупругости на основе трехпараметрического дифференциального операторного 
уравнения. 

Примем закон гармонического нагружения реологического тела в виде функции: 


&+2—&+471? /2= =4(), (1) 
т 


где 4 (+) — функция внешней нагрузки; = — деформация тела; ЕЛ] — коэффициент затухания 
(отношение модуля упругости к коэффициенту вязкости); /, — собственная частота и 
период незатухающих колебаний. 

Функцию внешнего воздействия в правой части уравнения (1) можно представить 
множеством способов, исходя из удобства рассмотрения задачи. В любом случае функция 4 (+) 
— это произведение: 

4(1)= КМ, (2) 
где № - параметр нагрузки; К - константы, характеризующие объект воздействия нагрузки. 

Наиболее удобными можно считать следующие способы выражения функции +(6), иС- 
пользуя основные положения ре И Е. 

















+) в == т т 9580; (3) 
Е 
О их ОО >= РО (4) 
1 - м а. 
„йо я О(® = ИИ 4$); (5) 
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+() = АР = АТ г АР (6) 
57 5Гр 

$ ( )==А7@= би АО - АМИ -А-АЙО; (7) 

о ра т 
4(#)= кат АТО = в” = И ь (8) 

1 

4(#)= Е. от 975 отд” (9) 
+()- 0-2 Б0=5Бо= об: (10) 
++, = ы „6(0= ро - ›5(0; (11) 


4(#)= ГИ = ыы. 590 (12) 


ЛЕт ЛЕТр ЛЕ р 
В соотношениях (3)-(12) введены следующие а Р - жесткость тела; 





72) — поглощенная доза излучения; Ь — МОЩНОСТЬ поглощенной дозы излучения, 0(!) — на- 
грузка, выраженная через напряжения; р(!) - нагрузка, выраженная через механический им- 
пульс; (0, О(®) — соответственно механическая энергия и тепловая; ЛТ(#) - изменение тем- 
пературы при нагреве тела; Р(!) — мощность механического воздействия; 5, Л5(г) — энтропия 
и изменение энтропии; С, с, ` Темплоемкости соответственно атомная и массовая, 
а, ).‚ — тепловое расширение и теплопроводность; 5,/,р, т — площадь, объем, плотность и 
масса тела; ^./, — эмпирическая константа и время релаксации; 5(1),5(#) — первая и вторая 
производные от действующей нагрузки; и — коэффициент вязкости. 
Функцию нагрузки (1) можно определить также с помощью операторного уравнения: 
(0) = рб рб + р›б =К(р)о (0. (13) 


Преобразуем общую форму записи операторного уравнения к виду, соответствующему 
уравнению (1), используя коэффициенты: 





р›6+2— р6+41? Ёрив =. (14) 
т 
Выразим искомые коэффициенты через вспомогательные величины: 
Р2=Ро, (15) 
Е 
В =2—Р, (16) 
т 
Ро =4п’ Лу ро. (17) 
Тогда: 
т 
= р» = : 18 
Р› = Р> Ея ( ) 
Е. ШБ 
=2— р =—, 19 
Р В Р А Ет (19) 
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* От 








= р 20 
ь 2^Е?т ро 
„ О 
Ро= 4п” о Ро Ее (21) 
Ет 
й р 
=——^_., 22 
Подставим полученные коэффициенты в (14): 
т. р. р 
6+ 6+— о=ч0. (23) 
ЛЕ^т ЛЕт Ет 


Полученное уравнение показывает, что приведения ее к наиболее компактному виду 
удобно использовать в качестве функции нагрузки первое выражение в (3). 
Выполним подстановку и преобразования уравнения (23): 





в АЕ (24) 
Ц Ц Ц 
ИЛИ 
6+46+/^20=426(0. (25) 


Выразим функцию нагружения через параметры системы, соответствующие левой час- 
ти уравнения (25), и уравнения для функции напряжений и ее производных, получим: 


б=В^?ое*. (26) 
Упростим выражение: 
зы 
(—е”^) = Веб , (27) 
0 
согласно которому определим: 
6(#) =Во,е”(1-е”'). (28) 


Полученное уравнение характеризует релаксацию нагрузки после ее установления. В 
случае, когда нагрузка изменяется по некоторому временному закону, полное выражение 
функции нагружения примет вид: 

$(0 =о(-Вое”(1-е”). (29) 

Такой вид функции нагружения наиболее компактен и удобен для анализа вязкоупру- 

гой среды. В большинстве случаев нагрузка (г) изменяется по гармоническому закону. Обо- 


значим максимальное значение возмущающей силы 4 


т! 


тогда, при Гармоническом законе 
изменения и в соответствии с условиями (3), можем записать: 


с з р 
+28 +008 =4, — ©0361 (30) 
Ет 
или в другом виде: 
а. 99 о 
&+2—&+ 0,5 =4, — с03 Е, (31) 
т № 


где ®=2л/ =2л/Т - угловая частота возмущающей силы и установившихся колебаний систе- 


МЫ, ть — частота и период нагружения. 


Решением уравнений вида (30)-(31) является уравнение [2]: 
==, с03(61 —0@), (32) 
где © — фазовый сдвиг резонатора (реологического тела) относительно возмущающей силы, 
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я © 206 204. 
@ = агсЁе 5 =авЮ = аг———— : (33) 
т\с —© © -© а 
т 


Таким образом, если нагрузка выражается через напряжения тела, то запаздывание 
колебаний составляет определенную выше величину угла механических потерь ЛАф, причем в 


большой степени зависящую от частоты вынуждающей силы © и коэффициента затухания 


р Е 1 
В=б=л=—=—, 1, >0. 
т 
Следовательно, в данном случае а =Аф . Тогда [3]: 


п” * 


[®] 





@ = Дф = асе = агс®—°. (34) 
Е 50 
Преобразуя уравнения (33) и (28), получаем: 
2ФЕ Е № Ц м 
Дф = агсе = агсе —— 1-е” ] = агс® — 1-е”). (35) 
пе уни ет) 


Согласно теории колебаний [2], амплитуда колебаний реологического тела изменяется 


от частоты возмущающей силы « и коэффициента вязкого трения В: 
1 


а 
7„(6)=Ч„ [и (8 -в*) бо? , (36) 
ИЛИ 
2] 2 
д. Ш 2 Е 
7„(6)=—” (%-®”) +4 Е) (37) 
т Ц 
Запишем это выражение с учетом эмпирического коэффициента: 
1 
ь 1 
ба (53 -о) ль? = (38) 
т 


Анализ уравнений (36), (37) показывает, что при ®-— ®, и при В О амплитуда мо- 


жет превысить допустимые для реологического тела значения, определяющие прочность тела. 
Динамические условия разрушения можно получить из условий прочности. При заданной воз- 


мущающей силе (1) и коэффициенте вязкого трения В амплитуда У, (6) является функцией 
только угловой частоты нагрузки. При ®«=®, достигается резонанс. Для определения резо- 
нансной частоты необходимо найти максимум функций (36)-(38) и приравнять первую произ- 


водную нулю, тогда 
62 
с = |8 5 =\/ 2 — 26? р (39) 
т 


Используя эмпирические константы и выражения связи с вязкоупругими характеристи- 
ками полимерных материалов, получаем: 





(40) 
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Параметром вязкого трения является коэффициент затухания бд. При малых значениях 
б величина У, (6) резко возрастает. Уравнения (36)-(38) справедливы и для случая статиче- 


ского нагружения тела при условии Ч, (®=0). 
Чтобы найти выражение для резонансной амплитуды, подставим (40) в (36)-(38): 











1 
ее еб (41) 
©) В [4 ры Е Во 2тбо 
1 2 
о (42) 
7„ (о) (о ^?) В [2 ы 


2т/\ю 2тбох ЗтЁЕ® 
Величина угла механических потерь Аф характеризует основные свойства реологиче- 


ского тела, материала или системы, а ее выражения (34)-(35) позволяют получить обобщен- 
ную связь между характеристиками тела и параметрами нагрузки. 

Выводы. Полученная динамическая модель вязкоупругого материала, конструкции или систе- 
мы при управляющем воздействии нагрузки, действующей по заданному закону в течение 
некоторого времени, позволяет осуществить эффективный анализ поведения реологического 
тела в целях оценки изменения его напряженно-деформированного состояния. 
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